Anomalous Behavior of the Contact Process with Aging 
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The effect of power-law aging on a contact process is studied by simulation and using a mean-field 
approach. We find that the system may approach its stationary state in a nontrivial, nonmonotonous 
way. For the particular value of the aging exponent, a = 1, we observe a rich set of behaviors: 
depending on the process parameters, the relaxation to the stationary state proceeds as 1/ \nt or 
via a power law with a nonuniversal exponent. Simulation results suggest that for < a < 1, the 
absorbing-state phase transition is in the universality class of directed percolation. 
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The study of systems with absorbing-state phase tran- 
sitions has been a topic of intense research in recent years 
Among these systems, the contact process (CP), 
introduced by T.E. Harris Q as a model for spreading of 
disease, has become the prototype model ||-^. In fact, it 
is a dynamical version of directed percolation . 

Recently, much effort has been directed toward general- 
izing the original models of absorbing-state phase tran- 
sitions by introducing many absorbing states, and to- 
ward finding paths to self-organized critical phenomena 
|p7|-Jl9[. Therefore, one can hardly expect that an ordi- 
nary CP will reveal new striking feartures. Nevertheless, 
in the present Letter, we demonstrate that even the sim- 
ple CP with a single absorbing state shows anomalous 
behavior if one introduces aging into the model. One 
should note that effects of aging, in the particular case 
of spreading of diceases were considered by Bernoulli |2^] 
long time ago. One may find a great number of models 
of evolution dynamics with aging in the book by Hop- 
pensteadt |21|. 

We show that the system with power-law aging, be- 
fore it approaches its stationary state, behaves non- 
monotonously: the time derivative of the particle density 
may change sign during the evolution. The case of the ag- 
ing exponent, a, equal to one is especially interesting. In 
this marginal situation, in a mean-field approach, we find 
different types of behavior. Depending on the parame- 
ters of the system, the density of particles changes at long 
times according to a power law with nonuniversal expo- 
nents, or proportional to 1/lnt. Such slow relaxation is 
certainly not typical of absorbing-state phase transitions 
with a unique absorbing state. At their critical point, 
relaxation normally follows a power law j^J^. 

We simulate the aging contact process (ACP) using 
sequential dynamics. Each site of the ID lattice may 
be empty or filled by one particle. At each unit of time 
a particle is chosen at random. We decide to annihilate 
this particle (a) or create a new one (b) with probabilities 
Pais) = 1/[1 + 1)"] or P,{s) = p{s + 1)"/[1 +p(s + 
1)"] correspondingly. Here, s is the age of the particle 
and a is the aging exponent. In the case of annihilation, 
(a), the particle is deleted and we proceed to the next 



step increasing time by one unit. In case (b), we chose 
one of two nearest neighbour sites with equal probability. 
If this site is filled, creation is impossible, and we increase 
time and proceed to the next step. If it is empty, we cre- 
ate a new particle at this site, increase time, and pass to 
the next step. Thus, as compared with the ordinary CP, 
we introduce the age-dependent probabilities Pais) and 
Pc(s). For a > 0, the case studied here, Pa{s) tends to 
zero and Pc{s) approaches 1 at long times. 

The process is started with a random configuration of 
particles. Simulation times are about 10* Monte-Carlo 
steps. The lattice size is taken up to 10^ sites, to mini- 
mize the fluctuations of the particle density, n{t), which 
is the main quantity of interest. The simulation results 
are presented in Figs. 1-3. First, one sees that at each 
a value studied, nit) may behave nonmonotonically, de- 
pending on the initial conditions and parameters of the 
process. If, for instance, the stationary state is ri(oo) — 1 
(see Figs. 2 and 3), the density may first decay nearly to 
zero, stay in this range a long time, and only approach 
the n(oo) = 1 much later. Frequently, it is only with 
long-time simulations that we can observe the final stage 
of the relaxation. 

Second, in the case of < a < 1, the critical charac- 
teristics are the same as for the contact process without 
aging, i.e., a = 0. In Fig. 1, we show the stationary parti- 
cle density as a function of the deviation from the critical 
point, Pc, for two different values of the aging exponent, 
a = 0.25 and 0.75. In both cases the data suggest that 
the critical exponent P is not a-dependent, and agrees 
with the directed percolation value, /3 ~ 0.277. When 
a > 1, the phase transition is absent and the only possi- 
ble stationary state is n(oo) — 1; the system exhibits an 
exponential approach to this limit (see. Fig. 2). 

In the marginal case, a — 1, we observe a complicated, 
slow evolution with several temporal regimes character- 
ized by different behaviors. Due to the slow relaxation, it 
is difficult to fix the final state. In Fig. 3 we present the 
results for the temporal evolution of the particle density 
for several values of the parameter p. Depending on p, 
different types of the evolution are observed. 

Let us present the simplest possible mean-field descrip- 



1 



tion of the ACP. While the model is even simpler than 
that used in our simulations, it is reasonable to expect 
that it can explain the principal features observed in the 
latter. As in the ordinary CP, each site of a lattice may 
be filled by one particle or be empty. Each particle has 
its age, s. At each increment of time, all lattice sites are 
updated according to the following rules, (i) The proba- 
bility for a particle of age s to die at the next instant is 
Ps, irrespective of its environment, (ii) The probability 
for a particle to survive until the next update is therefore 
1 ~ Ps. (iii) The probability that a particle will be born 
at an empty site in the next instant is (l/ns) 
where is the total number of nearest-neighbor sites, 
and n_F is the total number of particles at such sites. s{i) 
is the age of the particle at site i. (iv) The probability 
for an empty site to remain empty until the next up- 
date is therefore 1 — (I/'t-s) X]"=i 9s(i)- In this way we 
present a natural generalization of the CP, introducing 
age-dependent death and birth probabilities, Ps and g^, 
respectively. 

Let us derive the mean-field equations, introducing the 
following quantities: the total number of particles at time 
t, rit = J2\=o ^t,s, where at^s is the number of particles of 
the age s at time t. The initial condition is a 0,0 = '^(0), 
where n{0) is the initial number of particles. 

Since the only possibility for an old particle (s > 0) to 
survive is process (ii), the first equation is flf+i^^+i = 
(1 — ps)at^s, i-e., a fraction (1 — ps) of the particles 
with age s survive. One may verify the solution s — 

at-s,on«=o(l ~-P«) ' where n«=o = 1 by definition. 

We demonstrate the derivation of the second equation 
in the simplest case of only two nearest neigbouring sites 
(one may show that other coordination numbers lead to 
the same equation). Particles are only created via pro- 
cess (iii). Hence, 
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The previous equations describe completely the ACP in 
the mean- field approach. 

Applying X]l=o '^^^ ^^^^ mean- field equation one 
obtains rit+i-nt = at+ifi ~J2l=oP»"'t,s- Taking Eq. (||) 
into account, we immediately obtain the equation gener- 
alizing the usual mean-field equation for the CP without 
aging: 



Now we see that all quantities may be expressed 
through the at^Oi i-e., through the number of particles 
born at instant t. The equation for this variable follows 
fromEq. (|): 



at+1,0 



1 - ^at~s,o ]J(1 -Pu) 
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Formally speaking, Eq. (Q) together with the initial 
condition, define the solution of our problem. But the 
problem is still too difficult for analytical treatment. Let 
us pass to the continuum limit. Then the equations and 
the initial condition become 



7i(t) = / dsa{t,s) + a{t,t) , a(0,0) = n(0), (5) 
Jq 
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ds a{t, s)q{s) + a{t, t)q{t) 



/ ds a{t, s)p{s) + 
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a(t,0) = [l-n{t)] 
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One may check that the equations describing the contin- 
uous limit are self-consistent. Eq. (^) is the Von Foerster 
equation |22|, well known in populational dynamics. 
The solution of Eq. (||) is 



a{t, s) — a{t — s, 0) exp 



dup{u) 



(9) 



so a(t, i) = n(0) exp —J*dsp{s) . 

Here we consider only the simplest case oi q = const 
and an s-dependent p{s). From Eq. (^) we obtain the 
expression: a{t,0) = q[l — n{t)]n{t). This allows us to 
obtain, from Eq. (H), and accounting for Eq. (^), the 
following closed integral equation for the total number of 
particles: 



n{t) — n{0) exp 



dsp{s) 



(3) 



s=0 



q I c?s ?T,(s)[l — n(s)] exp 
/o 



dup{u) 



(10) 
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Differentiating Eq. (|T^) one obtains 
dn{t) 



at 



[q-p{t)]n{t)--qn'{t) 



qjds n{s) [1 — n{s)] [p{t) — p{t — s)] exp 





dup{u) 



(11) 



If we set p — const, Eq. (^l|) reduces to the usual equa- 
tion 1^,^ for the CP without aging. We write out the 
known results on linear relaxation for it for later compar- 
ison: ii p > q, n{oo) — and n(t) oc exp[— (p — q)t]] for 
p < q, ri(oo) = 1 —p/q and n{t) — n(oo) oc exp[— (q ~ p)t]; 
and at p — q, n{t) cc 1/t. 

The reasonable dependences for p{s) and q{s), which 
admit comparison with our simulations, are the follow- 
ing: p(s) decreases gradually as the particle age s in- 
creases, and q{s) increases with s increasing or is con- 
stant. Here, we consider the annihilation probability 
p{s) = c{s + io)~"j where a is the aging exponent. The 
constant c plays the same role as the parameter p^^ 
in our simulations. One can check that, in the partic- 
ular case of g = 0, i.e., when creation of particles is 
absent, Eq. is valid for any dimension, and ex- 

act. In this case, at < a < 1, n{t) approaches zero, 
n{t) = n(0) exp[iJ-"/(l - a)] exp[-(t + io)'-"/(l - a)], 
at a > 1, n{t) quickly approaches a constant value 
n(oo) = n(0) exp{— c^Q ^^/(a — 1)}, and, for a — 1, 
n{t) = 7i(0)ta/it + to). 

Eq. (M) may be easily solved by iteration. One may 
start, e.g., from n^^\t) — n(0). The resulting solutions 
n{t) are very similar to those obtained by numerical sim- 
ulation; we do not present these curves here, but only 
describe the results of the direct analysis of Eq. ( p^ ) . 

First of all, expanding Eq. (|lO|) for small t one finds 
that 



dn{0) 
dt 



-n(0){(7[l-n(0)]-rfo-"} . 



(12) 



(One may check that this relation is also valid for an s- 
dependent q{s) if one replaces q by q{0).) Hence, n{t) 
may increase or decrease from n(0) at short times de- 
pending on a particular relation between the constants 
of the problem and on the initial condition. In princi- 
ple, the sign of the derivative is independent of the value 
n{oo). That leads to the observed nonmonotonic behav- 
ior of n{t). 

The nonmonotonous behaviour observerd in some 
regimes can be intuitively understood regarding that in 
the very begining the particles are by necessity young, so 
most of them easily die off, but later the survivors turn 
to be nearly immortal, and the lattice will be filled up 
with old-timers. 

Let us now consider the behavior of n{t) for different 
values of a at long times. (Note that all the following 



asymptotes may be obtained from the Laplace transform 
of Eq. ( |lo| ) linearized near the corresponding stationary 
state.) 

Regime /, < a < 1. The situation is very similar to 
the CP without aging. There is a critical point c* {a, q, to) 
for each value of a and q. At c > c*(q;, tg), n(oo) = 
and n(t) has an exponcntional approach to the stationary 
state. At the critical point, c — c* (a, q,tQ), n(t) relaxes 
to zero by a power law, n{t) oc l/t. At c < c*(a, (7,io), 
at long times, n{t) approaches n(oo) > exponentially, 
n{t) — n{oo) oc exp{— g(a, c, q, to)t}. Here, g{a, c, to) is 
zero at the critical point, and also approaches zero when 
a approaches 1. 

Regime II, a > 1. For any c and g, n{t) approaches 
1 exponentially. The critical point (and scaling) are ab- 
sent. The reasons are the following. At a > 1, the kernel 
exp{— Jp duc{u + ^o) "} quickly decreases to a constant 
value, exp{— cip (a — 1)}, ast grows, and one can sub- 
stitute this constant into Eq. (p^). Taking the derivative 
of Eq. ( |lo|) and linearizing the resulting equation near 
n(oo) = 1 we get 

1 - n{t) oc exp |-(7exp[-cto ^""^V(a - 1)] i] ■ (13) 

Regime III, a = 1, the most intriguing case. Several 
types of critical behavior are realized for different values 
of c. 

(a) For c < 1 and any q, n{oo) = 1 and the long- 
time dependence is 1 — n{t) oc t^^^^'^\ Depending on 
the relation between c,q, and t^ [see Eq. (p^j, n{t) 
exhibits this dependence immediately or, on the con- 
trary, n{t) first decreases, remains near zero for some 
time, and then very slowly approaches 1. The kernel 
exp{— (isp(s)} = [{t -\- tQ)/tQ]~'^, so if we demand 
'n{t — > oo) ^ 1 in Eq. (p^), we immediately obtain such 
a form for the asymptote. 

(b) For c = 1 and any q, n{oo) — 1 and the long-time 
behavior is 1 — n{t) oc 1/lnt. 

(c) For 1 < c < 1 + qto, n(oo) = 1 — (c — l)/(gto) and 
•n{t) — n{oo) oc t^'^. 

(d) At c = 1 + qto, n{oo) = and n{t) oc l/t. The 
kernel is a quickly decreasing function as compared with 
ri(t), so the reason for such behavior is the same as for 
the l/t relaxation at the critical point at < a < 1. 

(e) At c > 1 -|- qto (note that the definition of p{s) 
imposes the additional restriction on possible values of c 
and c < tg), n{oo) = and n(t) oc t~'^. 

The exponents 1 — c and c are nonuniversal since c 
is simply a coefficient in the definition of the annihila- 
tion probability. Previously, nonuniversality was found 
in models with multiple absorbing states [p|JqJ^, [l8| , |T9| , p3t 
and in the process of spreading in media with memory 
[p^ . Nevertheless, the kind of nonuniversality observed 
here was not seen. 

One should note that the results for Regime III were 
obtained only in the frame of mean-field theory. We did 
not study the role of fluctuations, and our simulations 
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do not let us fix the final stationary state (except for the 
case n(oo) = 1), when the relaxation is slow, i.e., for 
a — 1. Note also that our mean- field approach is based 
on parallel dynamics while a sequentional one was used 
in the numerical simulations. This difference, however, 
seems to be not crucial for the CP with a single absorbing 
state. 

In summary, we have studied the influence of power- 
law aging on the contact process with a single absorbing 
state. Both simulation and mean-field theory reveal 
the nonmonotonic character of the relaxation. In the 
marginal case of a unit aging exponent, in the frame of 
mean-field theory, we found that the relaxation of the 
particle density to the stationary state proceeds by a 
power-law with a nonuniversal exponent, or even slower, 
i.e., proportional to 1/lnt. 
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FIG. 1. Log-log plot of the stationary density p vs. 
p — Pc for two different values of the aging exponent a, 0.25 
and 0.75. The critical points are pc(0.25) = 2.3419(3) and 
Pc(0.75) — 1.1014(3). The full line has a slope corresponding 
to the f3 exponent of DP in the 1+1 dimension. The inset 
shows the evolution of n{t) for different values of p (from top 
to below, p = 2.5, 2.4, 2.36, and 2.3423) and a = 0.25. 



FIG. 2. Evolution of n{t) for different values of p 
and initial densities, for a — 2. From top to below, 
p = 0.075, no = 0.01; p = 0.075, no = 0.25; p = 0.05, no = 0.05; 




FIG. 3. Temporal evolution of the particle density, n{t), 
for three different values of p and a — 1. From top to below, 
p = 0.65, 0.6667, and 0.75. 
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